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ESO-VLT and Spitzer spectroscopy of IRAS 05328-6827: a 
massive young stellar object in the Large Magellanic Cloud 
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ABSTRACT 

We present the first thermal-infrared spectra of an extra-galactic Young Stellar Object 
(YSO), IRAS 05328-6827 in the Hn region LHA120-N148 in the Large Magellanic 
Cloud. The observed and modelled spectral energy distribution reveals a massive YSO, 
M ~ 20 M Q , which is heavily-embedded and probably still accreting. The reduced dust 
content as a consequence of the lower metallicity of the LMC allows a unique view into 
this object, and together with a high C/O ratio may be responsible for the observed 
low abundance of water ice and relatively high abundances of methanol and CO2 ices. 

Key words: circumstellar matter - Stars: formation - Stars: pre-main sequence - 
ISM: individual objects: LHA 120-N 148 - Magellanic Clouds - infrared: stars 



1 INTRODUCTION 

Massive stars play an important role in galactic evolution. 
They deposit chemically enriched matter, momentum and 
ionizing photons into the interstellar medium (ISM), and 
their demise as supernovae may trigger or quench star for- 
mation in their vicinity. Because of their rapid evolution it 
is difficult to study the formation of massive stars. Until 
the protostar has become hot enough to create an (ultra- 
compact) Hn region it remains hidden from view inside its 
dust cocoon (e.g., Walmsley 1995; Garay & Lizano 1999). 
Studying the formation of massive stars in nearby, metal- 
poor environments has the advantage of reduced attenuation 
by dust, as the dust-to-gas ratio is a function of the metal 
abundance. An important scientific reason for studying mas- 
sive star formation at low metallicity is that it may be more 
representative of the early evolution of the Universe. 

Surveys for extra-galactic pre-main sequence stars have 
so far concentrated on prominent H 11 regions in the Magel- 
lanic Clouds. Optical and near-IR images reveal young stel- 
lar objects (YSOs) that have just emerged from their dust 
cocoon or are in the process of doing so (Rubio et al. 1998; 
Walborn et al. 1999; Brandner et al. 2001; Romaniello, Rob- 



berto & Panagia 2004). Embedded YSOs have been found 
by their dust emission at thermal-IR wavelengths or via 
maser emission at radio wavelengths (Whiteoak et al. 1983; 
Epchtein, Braz & Sevre 1984; van Loon & Zijlstra 2001; 
Maercker & Burton 2005). 

The mid-IR point source IRAS 05328-6827 is seen in 
the direction of the Large Magellanic Cloud (LMC), at 
RA= 5 h 32 m 38.59 s and Dec= -68°25'22.2" (J2000). On the 
basis of its near-IR colours and MSX 8.28 fim flux density 
it was selected as a target for our thermal-IR spectroscopic 
programme with the ESO-VLT and Spitzer Space Telescope 
of evolved stars in the Magellanic Clouds. Here we present 
evidence for it to be a massive YSO instead. We present and 
analyse the spectra and the spectral energy distribution, the 
first detailed investigation of an extra-galactic YSO. 



2 OBSERVATIONS 

2.1 Spitzer Space Telescope spectroscopy 

The InfraRed Spectrograph (IRS; Houck et al. 2004) on- 
board the Spitzer Space Telescope (Werner et al. 2004) was 
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used on 15 February 2005 to take a low-resolution spectrum 
(A/AA ~ 100) of IRAS 05328-6827 between 5.2 and ~ 36 
/xm, as part of Cycle 1 programme #3505. The star was ac- 
quired on the slit after peak-up in the red (22 /zm) on the 
target itself. The spectra were taken in staring mode, plac- 
ing the target on different positions on the slit to facilitate 
removal of sky emission. The spectral segments received a 
total integration of 28, 28, 60 and 420 seconds for respec- 
tively SL Order 2 (5.2-8.7 fim), SL Order 1 (7.4-14.5 /urn), 
LL Order 2 (14.0-21.3 (im) and LL Order 1 (19.5-38.0 /im). 
The spectral segments were combined by making corrections 
for pointing-induced throughput errors, normalizing to the 
best centered segment. The accuracy of the overall spec- 
tral shape is a few per cent. The reduction, extraction and 
calibration of the spectrum followed closely the procedures 
described in Sloan et al. (2004). 

2.2 ESO-VLT 3-4 (im spectroscopy 

The ISAAC instrument on the European Southern Obser- 
vatory (ESO) Very Large Telescope (VLT), Chile, was used 
on 7 December 2003 to obtain a long-slit spectrum between 
2.85 and 4.15 /im. The resolving power of A/AA ~ 600 was 
set by the ~ 0.6" seeing rather than the 2" slit width. 

The thermal-IR technique of chopping and nodding was 
used, with a throw of 10" and jittering by small amounts, to 
remove the background. The spectra were extracted using 
an optimal extraction algorithm. An internal Xe+Ar lamp 
was used for wavelength calibration The relative spectral 
response was calibrated by dividing by the spectrum of the 
B-type standard star HIP 020020. This removed most of the 
telluric absorption lines but introduced artificial emission 
features due to the photospheric lines of HIP 020020. The 
spectrum of IRAS 05328—6827 was therefore multiplied by 
a hot blackbody continuum with Gausian-shaped absorption 
lines of Bra 4.052, Pfy 3.741, Pf<5 3.297, Pfe 3.039 and PfC 
2.873 ^m scaled to match those in HIP 020020. 

2.3 Photometry 

HKL photometry of IRAS 05328-6827 obtained at the 
South African Astronomical Observatory (SAAO) was pub- 
lished by van Loon (2000). J HKL imaging photometry 
was obtained by us (PRW) at Mount Stromlo and Sid- 
ing Springs Observatory (MSSSO) using the ANU 2.3m 
telescope with CASPIR (McGregor 1994; Wood, Habing 
& McGregor 1998). JHK S photometry is available from 
the 2-Micron All-Sky Survey (2MASS; Cohen, Wheaton & 
Megeath 2003). Additional L'-band photometry is derived 
from our ISAAC spectroscopy acquisition images. 

An 8.28 /J,m flux density is taken from version 2.3 
of the Mid-course Space experiment (MSX) Point Source 
Catalogue (Egan & Price 1996; Egan, Van Dyk & Price 
2001). We collected scans from the IRAS data server 
(http:/ /www. astro. rug.nl/IRAS-Server/ I to measure 12 and 
25 ^m flux densities, using the Groningen Image Processing 
SYstem (GIPSY) software with the scanaid tool. 

The IRS red peak-up image provided photometry at 22 
^m (between 18.5 and 26 /im). The individual, shifted inte- 
grations were combined to remove bad pixels, and aperture 
photometry was obtained assuming a calibration of 1.16 data 
numbers /iJy -1 for a 4-pixel aperture radius. 



Table 1. Photometry of IRAS 05328-6827 (flux densities from 
MSX, IRAS & IRS are given in Jy). 
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Figure 1. MSX image at 8.28 (im showing the location in the 
LMC of the star forming region LHA 120-N 148. 

All the available photometry is listed in Table 1. No 
colour corrections were applied to the photometry. 



3 RESULTS 

3.1 The environment of IRAS 05328 6827 

IRAS 05328-6827 is seen against N 148 A, the north-eastern 
part of the H II region LHA 120-N 148 (Fig. 1). Compared to 
30 Doradus, located almost a degree to the south-east, and 
the H II region LHA 120-N 44, to the north-west, N 148 is not 
very bright but it does extend almost half a degree in diam- 
eter. Optical images are suggestive of the presence of a large 
dark cloud immediately to the north. N 148 A seems to be 
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Figure 2. A POSS2 red image (left), which includes Ha emission, shows the H II region LHA 120-N 148 A and a compact cluster embedded 
in it. The cluster stands out better on the POSS2 IR image (centre). IRAS 05328—6827 only becomes visible at wavelengths beyond ~ 1 
fim and is one of the brighter stars in the MSSSO i\-band image. The IRS red peak-up image (right) shows an intense source at 22 (im. 
Coincidcntally, the blue peak-up image at 16 fj,m includes the cluster, IRAS 05330— 6826, and a stellar mid-IR source, MSX LMC 1780. 



disjoint from the rest of N 148 (Fig. 2). There is no evidence 
for severe interstellar reddening of stars in the nebula. 

IRAS 05328—6827 is a bright object in the if-band and 
in the IRS red peak-up image at 22 /j,m (shown in Fig. 2 
in vertical alignment with, and at an identical spatial scale 
as the optical and if-band images). The blue peak-up im- 
age at 16 (shown in Fig. 2 in its true position relative 
to the red peak-up image) happened to catch the two other 
bright mid-IR sources in this region, IRAS 05330— 6826 and 
MSX LMC 1780 (both are easily visible in Fig. 1). The na- 
ture of MSX LMC 1780 is unclear, but IRAS 05330-6826 co- 
incides with an anonymous cluster of stars embedded within 
the Ha nebulosity. The mid-IR emission might arise from 
circumstellar material of cluster pre-main sequence stars or 
from diffuse dust that is heated up by the young stars. 

Neither 12 COi_ (Israel et al. 1993) nor radio emis- 
sion (Marx, Dickey & Mebold 1997) is detected from 
IRAS 05328-6827. The lack of radio emission (at a level 
of 6 mjy at 1.4 GHz or 3 mjy at 2.4 GHz) rules out the 
possibility that IRAS 05328-6827 is an ultra-compact Hn 
region (cf. Zijlstra 1990) unless it is optically thick at those 
wavelengths. Interestingly, very weak CO and radio emission 
at 1.4 GHz is detected from IRAS 05330-6826 . 



3.2 Modelling the spectral energy distribution 

The spectral energy distribution (SED) was reproduced with 
the dust radiative transfer model dusty (Ivezic, Nenkova & 
Elitzur 1999), as a first attempt at deriving quantitative 
information about this YSO. The fit is shown along with 
the observed SED in Fig. 3. 

We placed a blackbody with an effective temperature 
of 30,000 K in the centre of a spherically symmetric dust 
cloud. For a distance to the LMC of 50 kpc the integrated 
SED yields a bolometric luminosity of 40,000 L©. This would 
correspond to a late-O type, 20 Mq zero-age main-sequence 
star (Hanson, Howarth & Conti 1997). Although the model 
does not fit the SED in perfect detail, the bolometric lumi- 
nosity is recovered with about 10 per cent accuracy. 

For the radial profile of the density in the dust cloud we 



adopted a single power law of the form p(r) oc r~ 0,8 . Such 
a flat profile is consistent with an accreting YSO (Osorio, 
Lizano & D'Alessio 1999) or with a flattened geometry much 
alike that of a disc (e.g., Whitney et al. 2003). The dust 
cloud extends from an inner radius of a few hundred AU, 
at which the dust temperature reaches 800 K, to an outer 
radius of about one parsec. The model fit uses oxygen-rich 
cold silicates (Ossenkopf, Henning & Mathis 1992) with a 
uniform grain size of 0.07 /im, and requires an extinction of 
Ay = 15 mag. 

An improved fit of a spherically symmetric model to the 
SED cannot be achieved without invoking a two-component 
model, with a warm inner dust shell and a cold outer en- 
velope. However, the geometry in massive YSOs is almost 
certainly not spherically symmetric. The present data are 
unable to constrain detailed models. 



3.3 Spectral analysis: ionization, ices and dust 

The 2.9-4.1 ^m VLT/ISAAC spectrum and 5.2-36 
Spitzer/IRS spectrum are displayed in Fig. 4. We define ad 
hoc continua (dotted lines) and compute the optical depths 
with respect to this pseudo-continuum (Fig. 5). Although 
this has little meaning for the dust features it facilitates 
computing column densities for the molecular ices. 

The absence of any atomic emission lines (Fig. 4) makes 
it unlikely that an advanced ionized region is present. We 
estimate upper limits to individual hydrogen and helium 
emission line fluxes of Fii ne < 1 x 10~ 17 W m~ 2 , which 
corresponds to a luminosity of Lii ne < 3 x 10 26 W (=1.5 
Lq) or iphoton < 6 x 10 45 s _1 . In a Case B scenario, in 
which all lines except Lya are optically thin, the strongest 
line in the 2.9-4.1 /im region is expected to be Bra. The 
Bra:Br/y intensity ratio is ~ 3 (Storey & Hummer 1995) 
and by scaling we obtain an upper limit of 10 26 W to the 
Br7 luminosity. This is a typical value for detected Br7 in 
galactic massive YSOs but it is not unusual for it to be 
weaker (Ishii et al. 2001). 

The heavily smoothed VLT/ISAAC spectrum (Fig. 5) 
shows absorption by H2O (water) ice around 3.1 pm. There 
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Figure 3. Spectral energy distribution of IRAS 05328-6827 as Figure 4. VLT/ISAAC (top) and Spitzer/IRS (bottom) spectra 

observed via photometry and spectroscopy between 1 and 36 fim, of IRAS 05328-6827. Ad hoc continua (dotted) are drawn to aid 

and the DUSTY model fit (sec text). in analysing the absorption features (sec Fig. 5). 



is also evidence for absorption by CH 3 OH (methanol) ice at 

3.4 and 3.5 jj,m. The 3.4 jj,m absorption in galactic YSOs 
often peaks at 3.47 fim, but a peak around 3.40 /tm as seen 
in IRAS 05328-6827 is predicted if relatively little H 2 and 
H2CO (formaldehyde) ice is present (Schutte et al. 1996). 
Related features might be detected at 3.85 and 3.95 /jm 
(Dartois et al. 1999). The optical depth ratio at 3.1 and 3.47 
/jm is with ~ 3 an order of magnitude smaller than observed 
in galactic YSOs (Brooke, Sellgren & Geballe 1999). 

There is clear evidence for absorption at 15.2 pm due to 
CO2 ice (Fig. 5); both the position and shape are identical 
to that seen in the galactic massive YSO W 33 A (Boogert et 
al. 2004). CO2 ice is expected in a low-mass embedded YSO 
(Bergin et al. 2005), but it is often relatively weak in massive 
YSOs (Watson et al. 2004) unless significant processing has 
occurred (Gibb et al. 2004). 

The integrated optical depths J rdv (wavenumber v) 
of the ice features are converted into column densities 
by division through by the integrated absorbance values 
from Hudgins et al. (1993). We thus obtain iV(H 2 0) ~ 

1.5 x 10 17 cm~ 2 , and relative abundances of approximately 
H20:CH 3 OH:C02=2:3:l. Although these abundances come 
with considerable uncertainty, CH3OH seems to be more 
abundant than CO2 and H2O is clearly under-abundant. 
The galactic massive YSOs W33A and AFGL7009S also 
have an exceptionally high CH3OH abundance, but only 
marginally higher than the CO2 abundance and a few times 
below that of H 2 (Dartois et al. 1999). IRAS 05328-6827 
may thus present an extreme case of processing of the ice 
mantles on the dust grains (cf. Gibb et al. 2004) . 

We speculate that the low metallicity in the LMC may 
lead to a less tenuous dust environment surrounding YSOs, 
allowing interstellar photons to penetrate deeper and to be 
scattered around within the extended envelope. The photons 
can then affect the ice coatings on the grain surfaces. The 
lower oxygen abundance and possibly higher C/O ratio in 



the ISM of the LMC may also have an effect on the chemical 
balance in YSOs. 

The observed opacity ratio of the 18 and 10 fjxa silicate 
components exceeds any that could be reproduced using the 
optical constants available to us. This ratio is known to be 
high at the low temperatures encountered in the ISM and 
star forming regions, as compared to the warm circumstellar 
envelopes of evolved stars (Ossenkopf et al. 1992; Suh 1999). 
It goes beyond the present analysis to design a grain with 
matching optical properties. 



4 CONCLUSIONS 

We have discovered a massive Young Stellar Object, 
IRAS 05328-6827, in the LMC Hn region LHA 120-N 148. 
We present spectra at 2.9-4.1 fj,m obtained with the ESO- 
VLT and at 5.2-36 /im obtained with the Spitzer Space Tele- 
scope: the first thermal-IR spectra of an extra-galactic YSO. 
This opens the way to studying the early stages of massive 
star formation under unique environmental conditions (e.g., 
at low metallicity) not encountered in the Milky Way. 

IRAS 05328-6827 is a luminous (40,000 L©), massive 
(20 Mq) YSO embedded in an optically thick dust cloud. Its 
SED is reminiscent of Class I or flat-spectrum sources, and 
shows absorption features of silicate dust and ices. There is 
no sign yet of an emerging Hn region. 

At the lower metallicity in the LMC, compared to star 
forming regions in the Milky Way (by about a factor 4, 
Koornneef 1982), the lower dust content has two effects: (i) 
it becomes easier for the observer to detect the inner, accret- 
ing part of the cloud core, and (ii) interstellar radiation can 
penetrate deeper into the cloud and affect the ice mantles 
on the dust grains. Together with a higher carbomoxygen 
ratio this might be responsible for the observed low H2O ice 
abundance and high CH3OH abundance. 
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Figure 5. VLT/ISAAC (top, heavily smoothed) and Spitzcr/IRS 
(bottom) spectra of IRAS 05328—6827, expressed in optical depth 
with respect to the continuum drawn in Fig. 4. There is clear 
evidence for absorption by ices of H2O, CH3OH and CO2. 
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